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The vortex lattice structure in a d^2_j,2-wave superconductor is investigated near the upper 
critical magnetic field in the framework of the Ginzburg Landau theory extended by including the 
correction terms such as the higher order derivatives derived from the Gor'kov equation. On lowering 
, ^ . temperature, the unit cell shape of the vortex lattice gradually varies from a regular triangular lattice 

^ to a square lattice through the shape of an isosceles triangle. As for the orientation of the vortex 

, lattice, the base of an isosceles triangle is along the a axis or the b axis of the crystal. The fourfold 

^ ■ symmetric structure around a vortex core is also studied in the vortex lattice case. It is noted 

' that these characteristic features appear even in the case the induced s-wave order parameter is 

absent around the vortex of the d^2_j,2-wave superconductivity. We also investigate the effect of 
the induced s-wave order parameter. It enhances (suppresses) these characteristic features of the 
f-H I d^2_j,2-wave superconductor when the s-wave component of the interaction is attractive (repulsive). 
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Much attention has been focused on a vortex structure in a high-Tc superconductor. It is of great interest how 
the vortex structure of the high-Tc superconductor is different from that of a conventional superconductor. Various 
' investigations of the flux line lattice (or vortex lattice) for conventional superconductors are compiled in Ref. 
P5 including the gap and Fermi surface anisotropy effects. By a number of experimental and theoretical investigations, it 
O is concluded that the symmetry of this superconductor is most likely to be (i2,2_y2-wave. Therefore, one of the points is 
how the vortex structure of a da;2„y2-wave superconductor is different from that of an isotropic s-wave superconductor. 

The amplitude of the dx2-y2-wave pairing, i.e., — ky in momentum space, has fourfold symmetry for the rotation 
about the c-axis. We expect that, reflecting this symmetry, the core structure of an isolated single vortex may break 
the cylindrical symmetry and show fourfold symmetry in a dx^-y^-wsLYe superconductor, and the fourfold symmetry 
of the core structure may affect the vortex lattice structure. 

Keimer et alu observed an oblique lattice by a small-angle neutron scattering (SANS) study of the vortex lattice 
in YBa2Cu307 in a magnetic field region of 0.5T< H <5T applied parallel to the c axis. They reported that the 
vortex lattice has an angle of 73° between the two primitive vectors and is oriented such that the nearest-neighbor 
] direction of vortices makes an angle of 45° with the a axis. rXtie oblique lattice with 77° was also found by scanning 
tunneling microscopy (STM) at 6T by Maggio-Aprile et alu'U They observed the elliptic-shaped STM image of the 
4-^ vortex core, and concluded that this oblique lattice cannot be explained by considering only the effect of the injfcifinsic 

a in-planc anisotropy, that is, the difference of the cohereiice lengths between the a axis and the b axis directions .ala On 
the other hand, at low fields of 35 and 65G, Dolan et aln reported the regular triangular lattice by the Bitter-pattern 
. observation of the vortices, where the slight distortion of the vortex lattice is explained by the anisotropy of the 
' coherence lengths between the a axis and the b axis directions, ft is suggested that this deformation from the regular 
Q , triangular lattice in the high field region of a dx^-y^-wsLve superconductor is due to the effect of the fourfold symmetric 
O ' vortex core structure. 

^ I The fourfold symmetric vortex core structure in a c?2,2_j^2-wave superconductor was so,fai; derived theoretically when 
•"^ ■ the s-wave component is induced ai»und a vortex Q£.cLj!_j,2-wave superconductivity.LTtil This mixing scenario was 
" mainly put forth by Berlinsky et a/.Lla and Ren et al\12~i3 in the two-component Ginzburg Landau (GL) theory for s- 
^ . and d-wave superconductivity. According to the consideration based on the two-component GL theory, it is possible 
■ - - ' that the s-wave component is coupled with the d-wave component through the gradient terms. Therefore, the s-wave 
component may be induced when .thp d-wave order parameter spatially varies, such as near the vortex or interface 
under certain restricted conditions.E2JE3 The induced s-wave component around the vortex is fourfold symmetric. The 
resulting vortex structure in the d-wave superconductor, therefore, exhibits fourfold symmetry. 

However, it is noted that the amplitude of the induced s-wave order parameter strongly depends on the detail 
of the pairing interaction. In the weak coupling theory of the continuum model, for example, we define that the 
s-wave component of the order parameter is proportional to V^^, the s-wave component of the pairing interaction.Efl 
Therefore, when \ Vs \ is negligibly small compared with the dominant (i^2_j,2-wave component of the pairing interaction. 
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the induced s-wave order parameter can be neglected. On the other hand, it was reported that in the BogoHubov-de 
Gennes theory of the extended Hubbard model the apiplitude of the induced s-wave order parameter is supppessed 
when the on-site repulsion of the interaction increases.113 In the numerical calculation by Wang and MacDonald,ll3 the 
induced s-wave order parameter is small, on the order of a few percents of the dominant d^2_y2-wave component. 

For a moment, we consider the limit where the induced s-wave order parameter is negligibly small compared with 
the dominant d.j;2_y2-wave component. (We denote this limit as "pure (ij.2_j,2-wave" .) In this case, the two-component 
GL equations are reduced to a one-component GL equation for only a dx2-y2-wave component, which is the same 
form as that of the isotropic s-wave case within the conventional GL theory. Thus, the conventional GL theory gives 
the result that the vortex structure in a pure fi2,2_j,2-wave superconductor is exactly the same as that of an isotropic 
s-wave one, that is, the core-structure of a single vortex is circular symmetric and the vortex lattice is the triangular 
Abrikosov lattice structii¥G.U 

However, Ichioka et alxB showed that the core structure of a single vortex exhibits the fourfold symmetric structure 
even in a pure dj.2_y2-wave superconductor by a numerical calculation based on the quasi-classical Eilenberger theory, 
which can be applied even at low temperatures. In their calculation, the fourfold symmetric structure becomes clear 
on lowering temperature T. These facts indicate the following results. Strictly speaking, the conventional GL equation 
is valid only near the transition temperature Tc- Far from Tc, we have ixL.include several correction terms which is 
higher order in the small parameter ln(Tc/r) ~ 1 — Tc/T. Enomoto et altB showed that the correction terms consist 
of higher order terms of the order parameter and the higher derivatives (or non-local terms) which are neglected in the 
conventional GL theory. And they indeed succeeded in showing that the fourfold symmetric core structure of a single 
vortex in a pure dx2-y2-wave superconductor is reproduced by including these correction terms. This indicates that, 
in the study of the fourfold symmetric vortex core structure in a d^2_y2-wave superconductor, we have to consider the 
effect of the correction terms of the order \n{Tc/T) in addition to the effect of the induced s-wave order parameter. 
These correction terms may affect the structure of the vortex lattice. i-, |-, 

As for the effect of the induced s-wave order parameter on the vortex lattice structure, Berlinsky et alBu investi- 
gated it by using the two-component GL equations. As characteristic features of the dj.2_j,2-wave superconductor's 
vortex lattice, they suggested that the vortex lattice is deformed from a triangular lattice due to the induced s-wave 
order parameter and it leads to the double peak structure in the magnetic field distribution function P{h) as a func- 
tion of magnetic field h. In a pure (ia;2„j,2-wave superconductor, however, these features vanish within their theory 
since the induced s-wave order parameter is absent. On the other hand, the correction terms which introduce the 
fourfold symmetric core structure might give some effects on the vortex lattice structure even in the pure dj.2_y2-wave 
superconductor. In the study of the vortex lattice based on the GL framework, the effect of the correction terms in a 
dx2-y2-waYe superconductor has not been studied sufficiently. 

The purpose of this paper is to determine the orientation and the unit cell shape of the vortex lattice in the 
framework of the extended GL theory near the upper critical magnetic field Hc2- We also investigate the spatial 
variation of the vortex lattice structure, that is, the order parameter, the current and the induced magnetic field. 
In addition to the effect of the correction terms of the dx2-y2-wave superconductivity in the order hi{Tc/T), our 
calculation is performed by including also the effect of the induced s-wave component so that the contribution of both 
effects (the effects of the correction terms and of the induced s-wave component) can be estimated on the equal footing. 
Therefore, we reconstruct the two-component GL theory for s- and c?2,2_j^2-wave order parameters with including the 
correction terms of the order ln(Tc/T). The pure dr^2_y2--wave case is easily obtained in the limit where the amplitude 
of the s-wave order parameter is negligibly small. 

We consider the case of the vortex lattice under a magnetic field applied parallel to the c axis (or z axis) in the 
clean limit. Our consideration is restricted in the case of tetragonal symmetry described by the point group D4. 
Then, the a axis and the b axis of the crystal coordinate are equivalent each other. The orthorhombic symmetry case 
can be obtained by scaling the coordinate system, namely, the a axis direction and the b axis direction are made to 
be equivalent if the lengths are scaled by the coherence length of each direction. The Fermi surface is assumed to 
be two-dimensional, which is appropriate to high- Tc superconductors, and isotropic in order to clarify effects of the 
dj.2_y2-wave nature of the pair potential on the vortex lattice structure. The additional anisotropy coming from, e.g., 
the Fermi surface can be incorporated into our extended GL framework. 1— 11— 1 

The two-component GL equations were derived, for example, by Ren et aZJiil'Ej from the Gor'kov equation. However, 
their derivation is within the conventional GL framework. Here, to study the effect of correction terms of the order 
hi{Tc/T), we reconstruct the GL theory by including higher order terms of the order parameter and the higher 
derivatives other than terms in the conventional GL theory. In the pure dx2_y2-wave case, the present authors derived 
the GL equation with including these correction termSr&om the Gor'kov equation for the continuum model to study 
the fourfold symmetric core structure of a single vortex.113 In this paper, we derive the two-component GL equation for 
s- and dj.2_j,2-wave components from the Gor'kov equation with including the correction terms of the order hi{Tc/T). 
For simplicity of calculation, our derivation is performed in the weak coupling theory of the continuum model. As for 
other models, for example, Feder and Kallin derived the two-component GL equations from the extended Hubbard 
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model with including these higher order termsEJ In the GL equation derived from other model, the coefficient of each 
term is modified, but the derived terms have one-to-one correspondence to those of ours. Therefore, while our results 
may be modified quantitatively, our approach can be applied to the GL equation derived from other models and our 
results remain to be valid qualitatively. 

We calculate the dx2_y2--wave and the induced s-wave components of the order parameter near Hc2 by solving the 
linearized version of the derived GL equations. We also calculate the current and the induced magnetic field in the 
vortex lattice near Hc2 by using the current density expression which includes the contribution of the correction terms 
and the induced s-wave order parametjec. In our calculation of the vortex lattice structure, we use the magnetic Bloch 
state (or Eilenberger's basis function )EII as the basis function and expand various quantities in powers of \n{Tc/T). 
Then, we can determine the vortex lattice structure in terms of an analytically expressed form, which makes it easy 
to understand the detail of the properties of the vortex lattice, such as the dependence on the mixing of the s-wave 
component. 

As for the orientation and the unit cell shape of the vortex lattice, the equilibrium state is estimated by the minimum 
of the free energy. This was done by Abrikosov near Hc2 within the conventional GL theory for an isotropic s-wave 
superconductor .EZI Here, we have to extend his theory so that we can include the effects of the correction terms and of 
the induced s-wavgjCirder parameter. The extension about the effect of the induced s-wave order parameter was done 
by Berlinsky et alnu The extension about the effect of the correction terms wasjserformed by Takanaka in the case of 
isotropic s-wave pairing but anisotropic density of states at the Fermi surface£3 In this paper, we follow the method 
provided by Takanaka and extend his calculation to the case where the superconductivity is d^2_y2-wave symmetry 
and the induced s-wave order parametftr also exists. 

On the other hand. Won and MakiEja studied the vortex lattice structure near Hc2 numerically from the Gor'kov 
equation on the same stand point of ours. They suggested that the most stable configuration of the vortex lattice 
is a square lattice tilted by 7r/4 from the a axis at temperature T < 0.88Tc even in the pure dx2-y2-waMe case. In 
our interpretation, this deformation of the vortex lattice from a triangular lattice is due to the above-mentioned 
correction terms of the order ln(Tc/T), which are derived from the Gor'kov equation but neglected in the conventional 
GL theory. Our calculation is based on the GL expansion of the Gor'kov equation. Therefore, their numerical results 
are reproduced near Tc by our calculation as an analytically expressed form. We note that Won and Maki made 
some errors in their estimate of the free energy minimum. They estimated the free energy minimum by the Abrikosov 
parameter Pa which is originally used within the conventional GL theory. However, this expression of (3a itself should 
be modified by including the correction of the order ln(Tc/T). In their calculation, the free energy is estimated only 
in the case of a regular triangular lattice and a square lattice. In addition to them, the oblique lattice should be 
considered together. In our calculation, we derive the correct expression of the Abrikosov parameter. And, by using 
it, we estimate the free energy minimum by exhausting all the possible orientation and unit cell shape of the vortex 
lattice structure, including a shape of a scalene triangle. 

The rest of this paper is organized as follows. In Sec. ||, we construct the two-component GL equations with 
including the correction terms of the order ln{Tc/T). By using these equations, the da;2_y2-wave and the induced 
s-wave component of the order parameter are calculated near Hc2. In Sec. m, we calculate the current and the 
induced magnetic field in the vortex lattice state near Hc2 ■ The orientation and the unit cell shape of the equilibrium 
vortex lattice structure is determined by the free energy minimum in Sec. |l^. In Sec. 0, we present the spatial 
variation of the vortex lattice structure in the pure c?2,2_j,2-wave case. The summary and discussions are given in Sec. 
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II. ORDER PARAMETER NEAR Hc2 

We adopt the coordinate system shown in Fig. |l], where one of the unit vectors of the vortex lattice is along the x 
axis. Thus, the unit vectors of the vortex lattice are given by ri = {ax,0) and r2 = {(^ax,ay). The x axis forms an 
angle of 6*0 from the a axis of the crystal coordinate as shown in Fig. ^ Therefore the shape of the vortex lattice is 
characterized by ay/a^ and and the orientation of the vortex lattice is denoted by ^o- 

We consider the two-component GL theory for s- and (ia,2_y2-wave superconductivity. The order parameter (or pair 
potential) is given by 

A{r,0)^s{r)MO)+d{r)MO)- (2.1) 

Here, r — (x, y) is the center of mass coordinate of the Cooper pair. The direction of relative momentum of the 
Cooper pair is denoted by an angle 9 measured from the x axis in the ab plane. The symmetry function for s-wave 
{dx2 -y2-wa.ve) component (f>si0) i<Pdi(^)) is assumed to be real and to satisfy the normalization condition (0^) = 1 
ii'f'd} — Here, we write /(• • ■)d9/2TT = (• • •). In our presentation of figures, we use forms 
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</.,(0) = l, = \/2cos2(0-0o) 



(2.2) 



as an example. On the other hand, the pairing interaction is assumed to be separable as follows, 

v{e\e) = Vs4>s{o')4>s{e) + Vd4>d{s')M0), 



(2.3) 



where attractive interaction is treated as positive. In Eq. ( ^.3| ), the dx^-y^-'^SLve component Vd is dominant and the 
s-wave component Vg is smaller than Vd- 

The GL equation is derived as follows from the Gor'kov equation in the weak coupling theory of the continuum 
model (As for the detail of the derivation, see Appendix A of Ref. Now, ri(r)(f>(\i) and y(/)*(k)0(k') are denoted 
by A{r,e) and V{6,e'), respectively.), 

Np^A{r, e') = A„{V{e\ 9)A{r, 9)) + A^v^Vie', 0)(v • <i)^A{r, 6)) + A2{V{e', 0)|A(r, 0)pA(r, 9)) 
+ AivUV{9', 0)(v • q)*A(r, 9)) + 6Ai{V{9', 0)|A(r, 0)|%(r, 0)) 
+ 2A4vUV{9\ 0)[4{(v . q)2A(r, 0)}|A(r, 9)\' + A(r, 9f{{v ■ q)2A(r, 9)}* 

-2|(v • q)A(r, e)\^A{r, 9) + 3{(v • q)A(r, e)YA*{v, 9)]) + (higher order terms) (2.4) 

with V = v{9)/vp = (cos 9, sm9) and q = V/i + (27r/0o)A, where Np is the density of states at the Fermi surface, vp 
the Fermi velocity and 0o the flux quantum. The vector potential is given by A = (—Hqu, 0, 0) in the Landau gauge. 
When we consider the GL equation, the contribution of the internal magnetic field can be neglected in A compared 
with that of the external field Hq, since we consider the case of an extreme type II superconductor. In Eq. (^.4|), 
Aq = ln(2w£)7/7rr) with the Debye frequency ljd and the Euler constant 7, A2 = — /3/3 and A4 = a/3^/18, where 



62C(5) 
49C(3)2 



0.908. 



21C(3) 
167r2r2 



(2.5) 



with Riemann's ^-functions C(3) and C(5). The terms with A^ in Eq. (2.4) are the correction terms which are neglected 
within the conventional GL theory. 

\n{2uj-Dj / T^Tc) , we obtain 



Substituting Eqs. (|2T|) and (|J) into Eq. (|2J) and using the relation {NpVdy^ 
two-component GL equations for s- and (ia.2_j,2-wave components. 



Vo 



Vo 



2 



+ {fyl){2\d\h + dh*)} + 0(7') = 0, (2.6) 



-d+^{ (0^(v . q)2)d + (<^,0,(v . q)2)4 + A {(<^4^|d|2d + (<^2^2^(2|,|2d + ,2^*)} 
770 ?7o 3ryg 

-2d|(v • q)d|2 + 3{(v • ci)dYd*])\ + 0(7') = 0, 



(2.7) 



with 7 = ahi{Tc/T), the GL coherence length ^ — {/3i'|./61n(rc/r)}^/2 ^nd the energy gap in the absence of a 
magnetic field 770 — {ln(T'c/r)//3}^/2_ Because s(r) is in the order 0(7) as sh own later in Eq. ( 2.2(: ), the terms 
including s(r) are neglected among the correction terms of the order 0(7) in Eq. (2.7). The parameter Cs in Eq. (2.6) 
is defined by 

c-^ = a{{NpVs)-^ - {NpVd)-^ - HTc/T)} = a{{NpVs)-^ - {NpVd)-^} + 0(7), (2.8) 

which depends on the interaction parameters Vs and Vd- From Eq. (|2.8| ), > in the case Vs is attractive and 
T > Ts- Here, Ts defined by [NpVs)~^ — \'n-{2ujY)^ / ttTs) is the onset temperature of s-wave superconductivity in the 
hypothetical super-cooling state where dj.2„j,2-wave superconductivity does not occur even at T < Tc. Our study is 
restricted in the temperature region < T < Tc- On the other hand, < in the case Vs is repulsive. 

We note that, when the s-wave component Vs can be neglected in the pairing interaction (that is, in the limit 
Vs 0), Cs — > 0. Then, the pure d2;2_j^2-wave case is obtainecL In the case Cs = , w e recognize that s-wave 
component of the order parameter vanishes (s(r)=0) from Eq. ( ^.6| ). Therefore, Eq. (2.7) is reduced to the one- 
component GL equation for only the dx2-y2-wave component. There, the GL equation does not include the s-wave 
component, but includes the correction terms of the order 0(7). 
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In order to calculate the upper critical magnetic field Hc2 and the order parameter at Hc2 , we consider the linearized 
GL equations. The annihilation and the creation operators are, respectively, introduced as a = —e~^/'^{qx — iqy) and 
= ~e^^^^{qx + iqy), where e — inH/ipQ = Air/axCLy. Then, 



V • q 



(2.9) 



with v± = ± iVy)/'2 — By substituting Eq. (2_^) into Eqs. (2^) and (2/7), the linearized GL equations 

are written as 



1 2e£2 

—s + jcs—^{{(l)lv+v-)[aa%s+((t)s(l)dvl)a^d+{<f),(j)dv^_)a''^d}+0{j'^) = 0, 



(2.10) 



1 2e^2 



where [a™a'''"] means to take all the permutations of the product, such as [aa^ 

In the investigation of the vortex lattice state, it is convenient to use the magnetic Bloch statej^ 



(2.11) 

aa^ + a^a or [a^a''^^] = a'^a^^ + 



Mr\ro) = (e"n!)-i/2ai"V(r|ro|i)|*=o, 



where 



V'(i"|i"oW = (^~^) X!^^P ^r^{^^ '■^ ^ Vo+pay - 2t)^ + <^ - ipayi^x + xq + ^(P^x^ - Wo^^ 



(2.12) 



(2.13) 



Here, n is the index of the Landau level, and there are relations aipn =ry(pjl^n-i and ahpn^ \/ri^Vipn+i- The function 
'0o(r|ro) is the well-known Abrikosov solution of the vortex lattice. tIrEil In Eqs. ( ^.12 ) and (2.13), Tq — {xq^i/q) is 
related to the quasi- momentum of the Bloch state, and determines the position of the vortex center. In our presentation 
of figures, we set Tq = — (ri -f r2)/2 so that one of the vortex centers locates at the origin of the coordinate. 
The order parameters can be expressed as the linear combination of the magnetic Bloch states as follows, 



d{r) = ^ d„?/>„(r|ro), s(r) = ^ s„Vn(r|ro). 



(2.14) 



In the limit 7 ^ (T ^ Tc), Eqs. ( 2.10 ) and ( 2.11 ) are reduced to the form of the conventional GL equation. 
Therefore, the order parameters are given by the Abrikosov solution d{r) = doipo{r\ro) and s (r) = at T Tc. On 
lowering temperature below T^, the corrections of the order 0(7) are added to the Abrikosov solution. In order to 
investigate the corrections, the factors d„ and s„ (except for do) are expanded in powers of 7, 

d„=4«)-f7dW+0(72), (n^O) s„ = 4"^+74''+0(7'). (2.15) 

Now, the upper critical field is given by e = 47ri?c2/0o as the function of T. We also expand e in powers of 7, 

ee^^e(")e'+7e'''e' + 0(7'). (2.16) 

Equa tions ( |2.14| )-( [2l^ ) are substituted into Eqs. ( ^.1C| ) and ( |2.11| ). At the order of 0(1), we obtain Sn°' — from Eq. 
(^■lOD, and 



40) _ 2e("k'{{^li+v-){2n + 1)4") -|- (0.0^*^) x/(n + 2)(n + l)si°|2 + (</>.0d^!) \/n(n - l)^!"!^} 



from Eq. ( p. 111 ). Equation ( |2.17| ) gives 



(2.17) 



(2.18) 
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for n = 0, and dn'' = for n > 0. On the other hand, at the order of 0(7), we obtain 

siP + c,2e(")e' { {(bscl^dvl) Vin + 2)(n + + {cP,cPdvl)Vn{n - 1)4^1^] - 0, 



from Eq. ( 2.10 ), and 

{1 - 2e^"k'{c^li+i-){2n + 1)} 4^ - 2e^'k' {c^lv^i-)d'^^ 
-2e(")e' { {4>sC^dvl)V{n + 2) (n + l)^^ ^ + {4,,cf>dv'_)^n{n - l)si%} 

+2(e^'^fC^{{cf,lvl)^{n + 4)(n + 3){n + 2)(n + l)d^°l^ + {cl>lvt)Vn{n-l){n-2){n-3)d 



,(0) 



from Eq. (2.11). Equation (2.19) gives 



for n = 2, and sJr = forn 7^ 2. From Eq. ( |2.20D , we obtain 

2c 

for n = 0, 



2c, 



^2 2 



do 



for n = 4, and dn ^ = for n 7^ 0, 4. 
As a resuh, i7c2 is given as 



c2 



1 + 7^ + 0(7^ 



^0 2(</.2{)+f}_)e2 

with the use of e(i)/e('') in Eq. ( ^.22|) . The order parameters at Hc2 are given by 

d{r) = doVo(r|ro) + 74'V4(r|ro) + 0(7'), 



s(r)=74'V2(r|ro) + 0(7'), 



with the use of in Eq. (|2.23|) and s^^^ in Eq. (|2.2lD 



Here, the 0- integral is performed with adopting the forms of (pd and (ps in Eq. (2^). The values of 
in this paper are evaluated as follows. 



t)s4>dv\) = 



2-2 -2 • 
V,V_ 



1 



4V2 



1 




16' 






1 


dvl) 


^ 4V2' 



1 

— f 

32 



-Ai9a 



Therefore, Eqs. ( p.24D -(|2.26|) are reduced to 



e=^ = ^|l+7(2+c.)+0(7^) 



d{v) = doi Vo(r|ro) + ^7 Q + c, ) e-4.9«7^4(r|ro) + 0(7') [>, 



s(r) = do {-70,6-2'^° 7/^2 (r|ro) ^ o(^2)| 



(2.30) 



In the correction terms of the order 0(7) m Eqs. ( 2.28|)- ( 2.30| ), the terms without Cg comes from the correction 



terms of the GL equation, that is, the A4-terms in Eq. (2.4), which has been neglected so far in the conve ntiona l GL 
theory. The terms proportional to Cg comes from the induced s-wave order parameter. As shown in Eq. ( 2.30| ), the 
effect of the induced s-wave component appears in the or der (7). It corresponds to the result cZ(r)/s(r) l — T/T^ 



suggested by the two-component GL theory.! From Eq. ( 2.29| ), we recognize that the correction of the (i^2_j^2-wave 
superconductivity also appears in the same order of the induced s-wave component. 

For comparison, by exchanging (ps (pd and putting Cg = 0, we consider the case of an isotropic s-wave supercon- 
ductor without the s-d mixing. The corresponding ^-integrated values are evaluated as follows. 



]v+v-) 



1 

16' 



0, 



(2.31) 




{cs~Q) even in an isotropic s-wave 
the order parameter have the form 



In this case, from Eq. (2.31), i?c2 is written as the same form of Eq, 
superconductor within the order 0(7). On the other hand, from Eq. 
s(r) = soV'o(r|ro)i that is, the fourth Landau level function ip4^{v\vo) does not appear in an isotropic s-wave case. 
Then, it is noted that the appearance of V'4(r|ro) is the characteristic feature of the dx2 _y-2-WAYe superconductor's 



vortex lattice. From Eq. (2.29), we see that the induced s-wave order parameter enhances the weight of V'4(r|ro) term 
in the case Vs is attractive (c^ > 0). On the other hand, the weight of V'4(r|ro) term is suppressed in the case Vs is 
repulsive [cg < 0). The higher Landau level functions, tpsi'rlro), Vi2(r|ro), • • •, also appear in higher orders of 7. For 
example, the eighth Landau level function ips{r\ro) appears in the order 0(7^). 

As is seen from Eq. ( p. 28 ), Hc2 is enhanced by the induced s-wave order parameter in the case Vs is attractive, 
which is consistent with the result of Franz et al. (Fig. 9 of Ref. ||). On the other hand, Hc2 is suppressed in the case 
Vs is repulsive, which is consistent with the result of Won and Maki (Fig. 1 of Ref. p3[ ). 

fwe set ~ At with At = 1 — T/Tc in Eq. ( p.2g ), we obtain the result corresponding to that of Franz et 
There, the induced s-wave order parameter leads an upward curvature in the curve of Hc2{T) near Tc, since the 
factor of (Ai)^ is positive in Hc2- However, this upward curvature occurs even in the pure c?3.2_j,2-wave case (cs = 0) 
because of I7 in Eq. ( 2.2^ ), which is neglected in the study of Franz et al. but derived from the 0(7) correction of 
the dr^2_y2-wSiYe superconductivity. This upward curvature by I7 occurs even in an isotropic s-wave superconductor 
in this approximation. These unexpected results are modified as follows by a careful consideration. Exactly speaking, 
^-2 ^ r2ln(re/T) should be expanded as ^ At - |(Ai)2. Thus, we obtain 



c2 



At 



{Atr 



(2.32) 



Now, in the pure dj.2 _y2--w&ve case and the isotropic s-wave case (cs = 0), the factor of (At)^ is negative and the 
curve of Hc2{T) shows a downward curvature near T^. For Cg > |(1 — a), the factor of (At)^ changes to positive and 
the curve of Hc2{T) sho ws an upward curvature. 

As is seen from Eq. (2.21), the amplitude of the induced s-wave order parameter is determined by not only the 
parameter of the gradient coupling betwee n rf^ 2_„2- and s-wave order parameters, {(pgpdv^) : but also the parameter of 
the pairing interaction, Cg defined in Eq. ( |2.8| ). The amplitude of the induced s-wave order parameter is proportional 
to the produc t of b oth parameters. 

From Eq. (2.26), the spatial variation of the induced s-wave order parameter has the same structure as that of 
the second Landau level function V2(r|ro)- Figure § is the case of a square lattice {ay/ax = 0.5, C = 0.5). In Fig. 

^ (a), we show the amplitude of the induced s-wave order parameter, |s(r)/7S2^''| = |-!/'2(r|ro)|. The amplitude is 
suppressed near the vortex of the induced s-wave order parameter. To show it clearly, we schematically present the 




At the vortex center of the d,^ 



-wave 



position of the vortices and the winding number of their phase in Fig. 
superconductivity (its winding number is -1-1), the induced s-wave-jJi'dep-paijameter has a vortex with the winding 
number —1. It is consistent with the result of the single vortex caseE'aEyO'Ej'EJ Further, the s-wave order parameter 
has extra vortices with the winding number -f 2, which locates at the farthest points from the vortices of the dx2_y2- 
wave superconductivity. In total, we obtain the winding number -1-1 for the phase of the s-wave order parameter 
when we go round about the boundary of a unit cell in the vortex lattice. The amplitude and the winding number 
of the s-wave order parameter are not affected by the orientation 9o of the vortex lattice. Figure || is the case of an 
oblique lattice with 73° {ay /ax — 0.676, C = 0.5). In Fig. || (a), we show the amplitude of the induced s-wave order 
parameter, which is consistent with the result of Berlinsky et al. (Fig. 3 (b) in Ref. |^. The a axis and the b axis 
are exchanged each other in our calculation.) The corresponding position of the vortices and their winding number 
are shown in Fig. ^ (b). In the oblique lattice case, the vortex with the winding number +2 in the square lattice 
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case splits into two vortices with the winding number +1. In the limit of a triangular lattice, these vortices with the 
winding number +1 locates at the farthest points from the vortices of th e da;2_j,2-wave superconductivity. 
The spatial variation of the d2:2_j^2-wave order parameter in Eq. ( 2. 29] ) is presented later in Sec. 0. 



III. CURRENT AND MAGNETIC FIELD 



We calculate the current and its induced magnetic field around vortices near Hc2- Magnetic field H(r) is divided 
into an external field Hq = (0, 0, Hq) and an internal field h(r) which is induced by the current j(r), 



H(r) =Ho + h(r). 



(3.1) 



Amplitude of the order parameter is small near Hc2- Therefore, as for the current and magnetic field, we consider 
terms up to the order O(|dop). The expression for the current density is derived as follows from the Gor'kov equation 
with including the correction terms of the order 0(7) (As for the detail of the derivation, see Appendix A of Ref. p^.), 



j(r) = eV X H(r)/(0o/e') - x h(r)/(0o/e') 



27r r 



e(v{(v-q)A(r,0)}A*(r,0)) 



'7e'{(v{(v • q)3A(r, 0)}A*(r, 6)) + (v{(v • cO'A{r, 0)}{(v • q)A(r, 6)}*)] + 0{j' 



where terms with |A| (that is, terms in the order 0(|(io|'*)) are neglected. In Eq. (3.3) 



3c(/)o Stt^ 



where the GL parameter kql is given by 



77rC(3)A^; 



\e\) {UvpY 



With the substitution of Eqs. (pj]), (2.9) and (2.14) into Eq. (p^), the current density is written as 

jW = — 

ni,n2 

where 



(j)„i,„,V'ni(r|ro)V',t2(r|ro), 



7eC'{('^dVi)+)(v/("i + 3)(ni + 2)(ni + l)rf„i+3C + v/«2(n2 - l)(n2 - 2)d„,d;^_3 

+ V [ni + 2){ni + l)n2fi„i+2d*^-i + Vl*^! + 1)"2("2 - l)d„i+iC?,*2-2) 

+ (0^v{;3)(V"i('^i - - 2)d„,_3d;^ + V(»^2 + 3)(n2 + 2)(n2 + l)d„i<,+3 

+ \/ni{ni - l)(n2 + l)d„i_2rfj^2+i + V n-i("-2 + 2)(n2 + l)d„i_id*3_,_2) 

+ (?!)^vu2{)_)(^(3„^ + 2n2 + 4)V^I7TTd„i+i<^ + (2ni + 3^2 + l)^d„,d*„^_^ 

+ \/ini + 2)(ni + l)(n2 + l)d„j+2C?,*2+i + V"in2("-2 - l)fi„i-irf*2_2) 

+ {(l)lvv+v'i) (^{3ni + 2n2 + 1) V^Tdm-id^^ + (2ni + 3n2 + 4)\/?^2 + IdmCa+i 

+ v/ni(ni - l)n2d„i_2rf;2-i + + 1)K + 2)(n2 + l)d„i+id;2+2) } + ^(7^)- 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 
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Among the correction t erms of the order 0(7 ) in Eq. (3.7), we neglect the terms including s„, since s„ is in the order 
0(7) as shown in Eq. ( ^.26 ). Equation ( |3.7| ) satisfies a relation 



{(j) n 1,712 J 1 



(3.8) 



which ensures that j(r) is real . By su bstituting the results of e and the order parameters (i„ and s„ in Sec |l| [Eqs. 
(^.15D, (|2.16|), (|2.18D and (|2.2l|)-(|2.23D] into Eq. and by introducing 



we obtain 



1-27 



OiY)}\do\\ 



(3.9) 



(3.10) 



V6 



2c, 



2 \2 



0(7^) Mol 



(3.11) 



(i+)5,o = 7(0dW+w_)V54^^rfi5 



4 
30 



2c. 



{(l>s4>dV 



2 \2 



■0(7') Mol 



(3.12) 



(j+)2,i = 7 ((0s</'d*^)4'^rfS - (</'d«+«-) V2|do| 



V2 ( 



' 2 ^\{4>li+v 



2c, 



{(f)s(l>dvl){(l>s(t>dV-) 



blv+v^)^ 



Oh') Mop, 



(3.13) 



(j-)4,i = 7((/)rfW+w-)\/54^^rfo 



2c,, 



m2 \2 



+ 0(7') Mol^ 



The other components {j±)ni,n2 for ^1 ^ "-2 are zero with in the o rder 0(7). As a result, from Eqs. (3_^) and ( p.S| ), 
the current density is given as follows by using Eqs. (3.10) - (3.14), 



(3.14) 



= 2Re{(j+)i,oV'l'/'o + (i-)3,oV'3'i/'S + U+)5,0-ip5-ipo + {j+)2,l-ip2lp*i + {j-)4,li^4lpl}, 



(3.15) 



jy = 2Im{-(j+)i^o'01'0o + (i-)3,O'03'0o - (j+)5,oV'5V'o - ij+)2.l1p2lpl + (j- )44V'4'0t }• 



(3.16) 



The magnetic field induced by the current is obtained from the Maxwell equation. When we neglect terms of the 
order O(|(io|'*), the induced magnetic field h(r) has the following form. 



h(r) = J2 (h) V'«i(r|ro)C2('"|ro), 



(3.17) 



where (h)„j^„2 satisfies a relation 

(h) {(h)rii,n2}* 

for h(r) to be real. In our case, h(r) has only the ^-component, h(r) — (0,0, hz{r)). 



(3.18) 
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From the Maxwell equation (3.2), we have relations 



(j+)rai,ra2 = —-^\^~V^i^z)ni-l,n2 + V n2 + 1 (^Iz )ni ,ri2 + l | , 

(j-)ni,ri2 = — 2— + l(^2)ni + l,n2 ^ \/n2(/lz)«i ,n2-l | ■ 



We also have a relation 



Vni + l(j+)„i + l,„2 + V^(i-)"l-l,n2 - ^^^2 + l(j-)ni 



"2 + 1 



«2(j+)ni,n2-l = 



from divj = 0. Therefore, {hz)ni,n2 is written as follows from Eqs. (3.1C)-(3.14) and (3.19)-(3.21), 



2 



|^^^(j+)„,„_i=7(/.W)i,i + 0(f), 



where 



= - 



'"i^f-V op, «0) 



(^i^^)o,o 



(/^i^^)M 



(ft 



2(,^2^+i;_)2 
^/6 



{(f>s(f)dV+){<Ps4'dv'^ 



+ 2c, 



2 \2 



(3.19) 
(3.20) 
(3.21) 

(3.22) 

(3.23) 
(3.24) 

(3.25) 
(3.26) 
(3.27) 
(3.28) 



The other components {hz)ni,n2 for '^i ^ "-2 vanish within the order 0(7). As a result, the induced magnetic field is 
written as 



hAr) = {hfU^{hi'^)o.o} IVoP +7(/ii'^)i,i|V'i|' + 27Re{(ft«)4,oV'4V'o} +0(7' 
By using the 6'-integrated values in Eq. ( ^.27 ), hz{r) in Eq. ( 3.29| ) is reduced to 



1 - ^7 ( 1 + J ^ l^oP - ^7(1 + Cs)Hi\' + ^ 7 (1 + 2c.) Re {e-4'^°V^4^* } + 0(7') 



(3.29) 



(3.30) 



The spatial variation of the magnetic field in Eq. (3.30) is presented later in Sec. 

When we consider the case of an isotropic s-wave superconductor without the s-d mixing {(jjs (j)d, Cg — 0), hy 
using Eqs. (2.31) we obtain 



h^r) = ftfj (^1 - ^7) IV'oP - ^7lV^iP + 0(7') 



(3.31) 



where {hi^^)^^ vanishes. From the comparison between Eqs. (3.3C) and (3.31), it is concluded that the appearance 
of the V'4'00 term in the magnetic f ield d istribution is the characteristic feature of the d2:2_y2-wave superconductor's 
vortex lattice. As is seen from Eq. (3.30), the contribution of the ip4'4>Q term is enhanced by the induced s-wave order 
parameter in the case Vs is attractive (c. > 0). On the other hand, its contribution is suppressed in the case Vg is 
repulsive (c. < 0). 
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IV. ORIENTATION AND UNIT CELL SHAPE OF THE VORTEX LATTICE 



To know the equilibrium state of the vortex lattice, we search the free energy minimum about the orientation and 
the unit cell shape of the vorteaJattice in the vicinity of Hc2- The free energy F is calculated by following the method 
of AbrikosovtZl and Takanaka.E^I As a result, the free energy of the vortex lattice is written as follows (the detailed 
derivation is given in Appendix ^ , 



2(B-H,2)-h^' 



with Fo = Nf^'^tIq \n{Tc/T) and the magnetic flux density 



B = H(r) =Ho + h(r). (4.2) 



Here, we write J (• • •)(ir/^^ = (• • •). In Eq. (4.1), D4^[dc, sj is given by 



with 



o{^i)uT4 - , - 2|V'iPl^oP) + 0(7^) (4.3) 

3 dor 3((hiv+v) V / 



l^ol^ + 47Re <^ ^V4V'olV'oP > + 0(7'). (4.4) 



When a magnetic field has only the z component as in the present case, Eq. (|4.1|) is reduced to 



-1 



73-772^2 - IX /^2^2 \ 1 + ) ' (4-5) 



where we use the notations 



n( I /- a \ - F)i[dc,Sc] 1 hi 

C{ay/a.XM = - (4-6) 



(4.7) 



By u sing Eq. ( 3.29 ) and the orthogonalized condition ^/;„j^?/;*^ = i5„j,„2, the quantities in Eq. (45) and h?. in Eq. 
(4.6) are given as follows, 



(/ii'%,0 , (hi'') 



,(0) 



1 + 7' 



,(0) 



0(7'), 



(4.8) 



1 + 27 



,(0) 



|V'o|4 + 27 



(fei^^)M 



|i^op|i^iP + 47Re. 



(fei'^)4.0 



V'4V'olV'o| 



Oh' 



(4.9) 



In EqjM), 2R^ih,/h''°Y'C{ay/a^X, 60) corresponds to (2k —1)(3a of the Abrikosov's theory.tZI In the limit 7^0 
Eq. (O) is reduced to his result. 



0] 



The quantities IV'ol^i '04'0ol''/'op and |i/'ip|'0op are calculated as follows from Eq. ( 2.12 ), 

1/2 



iV'ol 



exp < — vr— (m^ + n^) > cos(27rCmn) 



(4.10) 
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X exp < — TT— (m + n ) + i2TiC,mn > , 
ax 



(4.11) 



IV'iPl^ol 



(4.12) 



These quantities depend on the parameters ay /ax and C characterizing the unit cell shape of the vortex lattice. 
Therefore, the free energy (4.5) depends on the unit cell shape of the vortex lattice through the quantities iV'ol'', 
V'4'0ol'0op and |i/;ip|V'op- For the triangular lattice case {ay/ax = •\/3/2, C = 0.5), V'4V'ol'/'op = 0. On the other 



hand, the orientation of the vortex lattice affects the free energy (^ through in Eq. ( |2.23| ) and {h 
in Eq. ( 3.2^^ as is seen from Eq. (2.27). Since hz does not depend on the orientation and the unit cel l sh ape, the 
free energy (4.5) has its minimum at the vortex lattice structure where C{ay/ax, Ci^^o) defined in Eq. (4.6) takes a 
minimum value. 



From Eqs. (4.3), (4.9) and (4.12), C(ay/ax,C,SQ) is written as 



C{ay/ax, C,0q) 



IK 



7c, 



{(t)s(t>dV+){(f>s4>dV- 



2 V3 



2^2 



7Re 



2c, 



V'4V'SIV'0| 



0(7^). 



(4.13) 



When we use the 6'-integrated values in Eq. (2.27), Eq. (4.13) is reduced to 

1 



C{ay/axX,So) = 



l-27-^U-:^7(ll + 7c,) 



V6 
4 



-^(1-^)7(1 



2c. 



)Re {6-4^^^4^-01^012} +0(7'). 



(4.14) 



The orientation of the vortex lattice is easily determined from Eq. (4.14). With respect to Oq, the minimum of 

C{ay I ax^COo) occurs at 



" i {^^^ ('/'4V'eSlV'op) - 71"! 



(4.15) 



for Cs > —0.5, where arg(i/)4?/;Q ItAoP) means the phase of the complex function ?/)4i/)g l^/iQp. Here, k » 1 since we 
consider the case of an extreme type II superconductor. For c, < —0.5, where Vg is strong repulsive, the minimum 
occurs at 



1 



-arg 



(?/'4?/'(SlV'0 



(4.16) 



Then, the orientation 6*0 is rotated by 7r/4 from that of the case Cg > —0.5. For both cases of Eqs. ( p^ ) and ( p^ , 
the minimum of C{ay/ax^ Ci ^0) as a function of is given by 



C {ay/ax,0 = mine„{C(ay/a^,C,6'o)} 



l-27-^^l-:^7(ll + 7c,) 



x/6 



- V 1 1 - ^ ) 7|1 + 2cM^rM^\ + 0(7'). 



(4.17) 



We numerically estimate the minimum of C{ay/ax,C) with respect to ay /ax ajid Q. The term of iV-ol^ prefers a 
triangular lattice, which is the well-known result in the conventional GL theory.EZl Then, at T = Tc (i.e., 7 = 6), 
triangular lattice is realized as the equilibrium vortex lattice. On the other hand, the term of — |i/'4''/'olV'op| prefers 
a square lattice. On lowering temperature from Tc (i.e., increasing 7) along the curve of Hc2{T), the contribution of 
— |'!/'4'0olV'oPI term increases. Then, the vortex lattice is deformed from a triangular lattice to a square lattice. There 
are three processes for this deformation, which are enumerated as (a), (b) and (c) in the following. 
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(a) On lowering temperature from Tc, the ratio ay/ax gradually varies from \/3/2 = 0.866 (triangular lattice) to 
0.5 (square lattice) with preserving C — 0.5. It means that the vortices form a shape of an isosceles triangle with 



BO = BA as schematically presented in Fig. |J. Since the function ■04V'olV'op in Eq. ( 4.11 ) is real and negative in 
this case, the orientation is given by 0o = for > —0.5. (In the following, our calculation is restricted in the case 
Cs > —0.5.) Therefore, the base OA of the isosceles triangle is along the a axis direction. 

(b) On lowering temperature from T^, the parameters [tty/ax, C) gradually vary from (%/3/2,0.5) [triangular lattice] 
to (1, 0) [square lattice] with keeping the shape of an isosceles triangle with OA = OB. In this process, varies from 
30° to 45° so that the base AB of the isosceles triangle is along the b direction. 

(c) On lowering temperature from T^, the parameters {ay/ax,0 gradually vary from (73/2,0.5) to (1,1) [square 
lattice] with keeping the shape of an isosceles triangle with AB = AO. In this process, Oq varies from —30° to —45° 
so that the base BO of an isosceles triangle is along the b direction. 

The cases (b) and (c) are related each other by the relation C{ay/ax, 1 — C: ~^o) = G{ay/ax, Ci ^o)- As a result, the 
cases (a), (b) and (c) are equivalent and give the same orientation and unit cell shape. The vortex lattice is deformed 
from a triangular lattice to a square lattice with keeping the shape of an isosceles triangle on lowering temperature 
from Tc. As for the orientation of the vortex lattice, the base of the isosceles triangle is along the a axis or the b axis 
direction. Therefore, at the lower temperature region, the equilibrium state of the-vortex lattice is a square lattice 
tilted by 7r/4 from the a axis, which is consistent with the result of Won and MakiE3 

To see the temperature dependence of the unit cell shape, we consider the case (a) in the limit k ^ 1. In this limit. 



Eq. (4.17) is reduced to 



V6 



Ciay/ax, C) = (1 - 27)|^o|4 - ^711 + 2c, 1 1^4^-01^01^1 + 0(7'), (4.18) 

where it is enough to consider only the contribution of the condensation energy. The ratio ay/ax in the equilibrium 
vortex lattice structure is evaluated from the minimum of C{ay/ax, C = 0-5) in Eq. ( 1.18| ). In Fig. ||, we present the 
temperature dependence of the ratio ay /ax along the curve of Hc2{T) for various mixing of the s-wave component, 
Cs = —0.2, —0.1, 0, 0.2, 0.4. As shown in Fig. H, the deformation from a triangular lattice starts just below on 
lowering temperature. In the pure dx2^y2-vjaYe case (c, — 0), the square lattice is realized for T < 0.858rc. In the 
case Vs is attractive (cs > 0), the square lattice is realized from higher temperatures (T < 0.888Tc for Cs = 0.2, 



T < 0.908Tc for = 0.4), since the contribution of — |V'4'0olV'o|^| in Eq. (4. IS) is enhanced due to the induced s-wave 
order parameter. On the other hand, in the case Vg is repulsive (cs < 0), the temperature region of the square lattice 
is suppressed (T < 0.837rc for = -0.1, T < 0.807rc for Cs = -0.2). 

When we consider t he ca se of an isotropic s-wave superconductor without the s-d mixing (^^ ^ (pd^ Cs — in Eq. 



( [4.13D ), by using Eq. ( |2.3l| ) we obtain 

C(a,K,C,eo) = 1^(1-27)-^ (^l-y7)}WF + 0(7')- 



(4.19) 



In this isotropic s-wave case, the contribution of '!/'4V'o IV'oj^ i^ absent from C{ay/ax, C? ^o)- Then, the minimum of {ipol"^ 
determines the equilibrium vortex lattice structure. A triangular lattice is, therefore, realized in all the temperature 
region. It is consistent with the well-kuowii result in the conventional GL theory.EJ 

In the calculation of Won and MakiJ£3o they obtain the equilibrium state by estimating the Abrikosov parameter 
/3a — (|A(r, ^?)|'')/((|A(r, 6')|2))2. It seems to be an analogy of the conventional GL theory. Within the order 0(7), 



we obtain 
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Pa = jdlVMol' + 0(7^) = \^Po\^ - ^7(1 + 2cs)Re{e'^''°^iro\M^} + 0(7^). (4.20) 



where it is noted that 1 — 27 in the first term of Eq. ( [4.14[ ) is changed to 1. However, Pa itself s hould be modified by 



the correction of the order 0(7). In our calculation, we use the function C{ay /ax, Ci^^o) in Eq. (4.14) instead of /3 



It changes the result quantitatively. Dominant terms neglected in the estimate by Pa come from the gradient term 
of the A'^-order such as {(0^(v • q)^)(i}|(i|^ in the GL equation. If we consider the correction of the non-local term 
((/i^(v • q)^)(i in the calculation of the A-linear order, in the estimate of the free energy we have to include also the 
gradient terms of the A^-order as correction, which appear within the same order of 7. 

While Won and Maki considered only the cases of a triangular lattice and a square lattice, they obtain the result 
that the square lattice is stabilized at T < O.SSTc in the pure dj,2_y2-case. This temperature is not so different from 
ours. In their result, as repulsion of the s-wave component increases in the interaction, the temperature region of the 
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square latticejis suppressed (Fig. 4 of Ref. ^3|). This tendency of the repulsivja case is consistent with our result. Won 
and MakiES'O suggested that the experimental result of an oblique latticeoQ can be explained by the square lattice 
structure obtained here and by the scaling of the coordinate system due to the difference of the coherence lengths 
between the a axis and the b axis. 



V. SPATIAL VARIATION OF THE VORTEX LATTICE STRUCTURE IN THE PURE Djf2_y2-WAVE 

CASE 



As clarified in previous sections, the characteristic features of the d^2_y2-wave superconductor's vortex lattice is 
produced by the appearance of the fourth Landau level function '04(r|ro) in the vortex lattice structure. Due to this 
effect, the vortex lattice is deformed from a triangular lattice to a square lattice as shown in Sec. IV. In this section, we 
consider the effect of ip4{r\rQ) on the spatial variation of the vortex lattice structure in a d^^-yi-wave superconductor 
by using the results of Sees. || and III. Here, we consider the case of the pure d^2_y2-wa,ve superconductor (cg = 0). 
Then, in the following characteristic features of the vortex lattice, there is no contribution of the induced s-wave order 
parameter. 

First, we consider the spatial variation of the da;2_j,2-wave order parameter, which is given by Eq. ( ^.2S| ). The 
contour plot of the amplitude, |c?(r)/c?o|, is presented in Fig. ^, where the core region of the vortex is darkly shaded. 
We consider the case of a square lattice at a lower temperature [Figs. |6| (a) and (b)] and the case of an oblique lattice 
at an intermediate temperature [Fig. ^ (c)]. The parameter 7 is chosen by following the result of Fig. ^. 

As for the case of a square lattice {ay/ax=0.5, C=0.5), we use 7=0.2 (T/Tc=0.802). In addition to the case of an 
equilibrium orientation 9o = 0° in Fig. ^ (a), we also show the case of an unstable orientation 9o = 45° in Fig. ^ 
(b) to clarify the contribution of 'i/;4(r|ro). In both figures, the a axis and the b axis are along the horizontal and 
vertical directions. (In Fig. ^ (a), the x axis and the y axis are also along the horizontal and vertical directions. 
In Fig. ^ (b), the x axis and the y axis are rotated by 45° from the case of Fig. ^ (a).) Due to the property of a 
dj:2_j,2-wave superconductor, the amplitude of the order parameter is slightly suppressed along the directiO|iL.p£ the 
a axis and the b axis around a vortex, which has been exhibited in the calculation of the single vortex case.ESES On 
the other hand, due to the property of the vortex lattice, the amplitude of the order parameter is slightly suppressed 
along the directions pf. nearest-neighbor vortices, which has been exhibited in the calculation of the isotropic s-wave 
superconductor case.a Therefore, when the a axis and the b axis agree with the nearest neighbor directions as in the 
case of Fig. ^ (b), the suppression along the nearest neighbor directions is enhanced and the fourfold symmetric vortex 
core structure is emphasized. On the other hand, when the a axis and the b axis are away from the nearest-neighbor 
directions as in the case of Fig. ^ (a), the suppression of both effects cancels each other. Therefore, the vortex core 
structure is reduced to a cylindrically symmetric one. It seems that the gain of the condensation energy in the case 
of Fig. H (a) is more efficient than that of Fig. ^ (b). Then, the equilibrium vortex lattice structure prefers the 
configuration of Fig. ^ (a). 

Figure H (c) is the case of an oblique lattice with 73° {ay/ax=0.676, C=0.5) for the equilibrium orientation = 0, 
where we use 7=0.115 (T/Tc=0.881). In this figure, the a axis and the b axis (and the x axis and the y axis) are along 
the horizontal and vertical directions. Also in this oblique lattice case for 6*0 = 0, the vortex core has cylindrically 
symmetric structure since the a axis and the b axis are away from the nearest-neighbor directions. 

Next, we consider the spatial variation of the current a nd th e induced magnetic field around vortices in eac h cas e 
of Fig. 0. The current is calculated from Eqs. ( 3.15 ) and ( 3.16| ), where we use the ^^-integrated values of Eq. ( 2.27 ). 
Figure^shows the contour plot for the amplitude of the current, |j(r)|/(7r-\/e^/2K^77^). The corresponding magnetic 

field is calculated from Eq. ( |3.30 ). Figure |^ shows the contour plot of the magnetic field, hz(v)/\hf^ \ . The contour 
lines are also interpreted as the stream lines of the current. 

The calculation of the single vortex lattice case shows that |j(r)| and hz{r) are fourfold symmetric around a vortex 
core in a dj.2_y2-wave superconductor .F^ri There, |j(r)| has four peaks around the vortex core at the direction of 45° 
from the a axis or the b axis. Reflecting this current distribution, the magnetic field extends along the a axis and the 
b axis directions around the vortex core. This is because the field is strongly screened by the induced current along 
the 45° directions. These fourfold symmetric structures of |j(r)| and hz{r) are also recognized in our calculation of 
the vortex lattice case, which becomes clear with increasing 7 (lowering temperature). 

In the square lattice case at a lower temperature (7 = 0.2), we compare the case of the equilibrium orientation 
9o = 0° in Figs. ^ (a) and ^ (a) with that of the unstable orientation 6*0 = 45° in Figs. (b) and [s] (b). For the 
equilibrium case = 0°, the screening current flows in the narrow region around each vortex compared with the 
unstable case 6*0 = 45°. Then, for 6*0 = 0°, the enhancement of hzi^) at each vortex core localizes in the narrow 
region, and hz{r) is reduced to be a flat distribution at the region apart from vortex cores. 
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Figure ^ (c) and ^ (c) are the oblique lattice case at the intermediate temperature (7 = 0.115). There, fourfold 
symmetric structure of |j(r)| and h.z{r) are not seen so clearly due to the small 7 and also due to the broken fourfold 
symmetry of the vortex lattice structure. 

Another way to investigate the spatial variation of hz{r) is to consider the magnetic field distribution function 
defined as 

P{h)^ f 5{h-h,{v))dv/a.,ay, (5.1) 

J unit cell 

which is related to the resonance line shape in the /iSR and NMR experiments. In Fig. ^, we show P{h) as a function 

of h/\h'"z \ for the equilibrium state of the vortex lattice at various temperatures; 7 =0 (a), 0.06 (b), 0.115 (c), 0.15 (d) 
and 0.2 (e), which correspond toTjTc =1, 0.936, 0.881, 0.848 and 0.802, respectively. The case of a pure dx^^y2-wscve 



superconductor given by Eq. (3.30) is presented by solid lines. The parameters of the vortex lattice are chosen by 
following the result of Fig. |. Then, C = 0.5, 6*0 =0 and ay/a^ =0.866 (a), 0.799 (b), 0.676 (c), 0.5 (d), 0.5 (e). In the 
case of a triangular lattice [line (a)] and a square lattice [lines (d) and (e)], P{h) has single peak structure. On the 



other hand, in the case of an oblique lattice 
are consistent with the result of Franz et al 



es (b) and (c)], P{h) has double peaks structure. These structures 
while the oblique lattice is derived from the induced s-wave order 
parameter in their theory. As explained by them, the double peak structure occurs in an oblique lattice because the 
saddle points of hz{r) between two vortices are not equivalent in the ri direction and in the r2 direction. 



For comparison, we also show P{h) of the isotropic s-wave case given by Eq. (3.31), which is presented by dotted 
lines in Fig. ^. In this s-wave case, the equilibrium state is a triangular lattice {( = 0.5, ay /ax = 0.866). We compare 
the case of a square lattice in a pure dx^-y^-^&^ve superconductor with that of a triangular lattice in a isotropic s-wave 
superconductor. At 7 = 0.16 [line (d)], the peak of the square lattice case locates at the tight side of the peak of the 
triangular lattice case, which is consistent with the result of the conventional GL theory.E3 In the d^^ -yi-wscve case, 
with increasing 7 (decreasing T), the lower edge of P{h) approaches the peak position and the peak height increases 
as is seen from the line (e). This change of the distribution reflects the result of Fig. ^ (a), where hz{v) has uniform 
distribution outside the core region. In the isotropic s-wave case, on the other hand, the distribution of P{h) does 
not so vary on lowering temperature. m ^ 

Some of the above-mentioned features are similar to that of Berlinsky et al.pa where the characteristic features of 
the vortex lattice are derived from the contribution of the induced s-wave order parameter. However, it is important to 
notice that the characteristic features suggested by our calculation occur even in t he pure da;2_j,2-wave case where the 



induced s-wave order parameter is absent. Its origin is the A4-terms in Eq. (2^), which is in the order 0{hi{Tc/T)) 
in the dimensionless form of the GL equation and neglected in the conventional GL equation. The contribution 
of these correction terms increases with lowering temperature from Tc and leads to the fourfold symmetric vortex 
core structure. It also induces the ■(/'4(r|ro) term in the vortex lattice structure and leads to the above-mentioned 
characteristic features of the vortex lattice in a d2;2_j^2-wave superconductor. 

If the induced s-wave order parameter cannot be neglected, the contribution of the ?/'4(r|ro) term increases for > 
and decreases for Cg < 0, as shown in the result of Sees. || and [II. Therefore, the above-mentioned characteristic 



features, which are produced by the V-'4(r|ro) term, are enhanced in the case Vs is attractive (c^ > 0) and suppressed 
in the case Vs is repulsive (c^ < 0). 



VI. SUMMARY AND DISCUSSIONS 



We investigate the vortex lattice structure in a d^^-y^-waNe superconductor near Hc2 in the framework of the 
extended GL theory. We determine the orientation and the unit cell shape of the vortex lattice. As for the unit cell 
shape, vortices form a shape of an isosceles triangle as shown in Fig. ^. The vertical angle gradually varies from 
60° to 90° when temperature is lowered from Tc. Then, a square lattice is realized at the low temperature region 
T < 0.86Tc for the pure dx2-y2-wave case (cg = 0). As for the orientation, the base of an isosceles triangle is along 
the a axis (or the b axis). 

We also calculate the spatial variation of the vortex structure, such as the (i2,2_j^2-wave order parameter, the 
current and the induced magnetic field. As temperature is enough lowered, the fourfold symmetric core structure 
clearly appears in our calculation of the vortex lattice case. This structure is consistent with that of the single vortex 
case. The amplitude of the current has four peaks around each vortex in the 45° direction from the a axis. Then, 
the induced magnetic field extends along the a axis and the b axis direction. As for the dx2-y2-wa.ve order parameter, 
the amplitude around the vortex is slightly suppressed along the a axis and the b axis direction due to the nature of 
the (ia;2_y2-wave superconductivity. On the other hand, the amplitude is slightly suppressed in the direction of the 
nearest- neighbor vortex. When the a axis and the b axis are away from the nearest-neighbor directions (i.e., the case 
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of the orientation 9o = 0°), the suppression of both effects cancels each other. Therefore, the vortex core structure is 
reduced to a cyhndrically symmetric one. This structure seems to most effectively gain the condensation energy. It 
is the reason why the equilibrium state prefers the orientation Oq = 0° . _ _ 

Some of the above-mentioned characteristic features of the vortex lattice are similar to that of Berlinsky et al.fTu 
where the origin of the characteristic features in their result is the mixing of the induced s-wave order parameter. 
However, it is noted that the characteristic features obtained in this paper appear even in the pure rfj.2_j,2-wave case 
where the induced s-wave order parameter is absent. The origin of these characteristic features is the higher order 
terms of the order parameter and the higher derivatives (or non-local terms) which are neglected in the conventional GL 
theory. These correction terms derived from the Gor'kov equation are higher order in the small parameter ln(rc/r). 
Therefore, their contribution vanishes at Tc, but becomes increasingly important upon lowering temperature from Tc- 
Some of these correction terms induce the contribution of the fourth Landau level function ■!/)4(r|ro) to the vortex 
lattice structure within the order 0(ln(Tc/T)). This contribution of V-'4(r|i"o) leads to the characteristic features of the 
vortex lattice in the d^2_y2-wave superconductor, such as the deformation from a triangular lattice and the fourfold 
symmetric structure around a vortex core. 

We also estimate the effect of the induced s-wave order parameter. The amplitude of the s-wave order par ameter 
depends on Vs, the s-wave component of the pairing interaction, through the parameter Cs defined in Eq. (2.8) in our 
model. When the contribution of the induced s-wave order parameter cannot be neglected, the factor of the ■(/'4(r|ro) 
term increases (decreases) for Cg > (cs < 0) in the expression of the vortex lattice structure. Then, the characteristic 
features, which is produced by the ?/'4(r|ro) term, are enhanced (suppressed) in the case Vs is attractive (repulsive). 

In the following, we give some comments. 

At further low temperatures, we have to consider further higher order derivatives in the GL equation. Then, the 
contribution of higher Landau level functions ■(/'8(r|ro), ■(/'i2(r|ro), • • •, appear in the vortex lattice structure, which 
are in the order 0([ln(Tc/T)]") (n > 2). The investigation of these higher order correction is now in progress. 

The fourfold symmetric behavior is also observed when magnetic field is applied in the direction parallel to the ab 
plane. When the direction of the field is rotated within the plane, the Hc2 (Refs. ^ and 28) and the torque (Ref. 
|29| ) experiments show the fourfold symmetric behavior. They are explai ned hy , the contribution of the higher order 
derivative terms (or non-local terms) in the GL equation as in this paper .cSe3H Also in this parallel field case, if the 
induced s-wa,ve oder parameter exists, the fourfold symmetric behavior is enhanced (suppressed) when Vs is attractive 
(repulsive) S3 

The amplitude of the induced s-wave order parameter depends on V,,, the s-wave component of the interaction. 



through the parameter Cg in Eq. (2.8). Then, next, we want to know the value of Vg appropriate for the high- Tc 
materials. However, it seems to be a difficult problem, which is deeply related to the pairing mechanism of the high- Tc 
superconductors. If the s-wave order parameter is defined as the on-site pairing, Vs may be strongly repulsive since 
the high- Tc materials are strongly correlated electron systems. However, s-wave order parameter is defined as the 
extended s-wave, we cannot exclude the possibility of the attractive Vs- The extended s-wave is derived from the 
non-local pairing between nearest neighbor atomic site, which is the same origin of the c?a;2_y2-wave pairing. 

We note that it is not clear whether the two-component GL theory can be applied to the case of strongly repulsive 
Vs- In the limit of strong repulsive Vs (c^ — > — oo), th e ind uced s-wave order parameter completely destroys the 
(ia;2_j^2-wave superconductivity, which is seen from Eq. ( 2.32 ) for example. In this strong repulsive case, we have to 
reconsider from the outset whether we can introduce the s-wave order parameter or not. 

So far, we consider the case of the isotropic Fermi surface. If the Fermi surface is anisotropic, the average (• • •) 
is modified to {■■■p{9)), where p{9) is the ^-dependent factor of the density of states on the Fermi surface and 
/ p{9)d9/2iT = 1. And |v| has the 6'-dependence. 

The origin of the fourfold symmetric core structure and the deformation from a triangular lattice is the appearance 
of V'4(r|ro) in the vortex lattice structure. The contribution of ■!/'4(r|ro) is proportional to the quantity in the 

pure (ia;2_y2-wave case. Even in the s-wave superconductor {(pd —f(j)s), when the structure of the Fermi surface or the 
energy gap {\(j)s\) have the anisotropic component with fourfold symmetry, the contribution of 1/14 (r|ro) appears since 
^ 0. Then, the fourfold symmetric core structure and the deformation from a triaugular lattice occur at lower 
temperature and higher magnetic field. The square lattice observed bji-.De Wilde et aZH3 in the STM image of the 
vortex lattice on LuNi2B2C seems to be this case. And Eskildsen et alS3 reported the deformation from a triangular 
lattice to a square lattice with increasing magnetic field by the SANS and magnetic decoration observation of the 
flux line lattice on ErNi2B2C, where the orientation and the unit cell shape of the vortex lattice in the deforming 
process agree with our results. In the dj:2_y2-wa,ve superconductor, \ {(j>'^v'tp)\ takes a large value since 0^ has fourfold 
symmetric node structure. In the s-wave superconductor, to introduce a large value of \{(j>'^v'^p)\ comparable to the 
dx2_y2-wave case, (j^Hv- j^p should have a large anisotropy with fourfold symmetry such as the form of (1 -|- cos40)/2. 



16 



ACKNOWLEDGMENTS 



We are grateful to K. Takanaka for informing us the detail of his calculation. We also would like to thank N. 
Hayashi and T. Sugiyama for fruitful discussions. One of the authors (M. I.) would like to thank the Japan Society 
for the Promotion of Science for financial support. 



APPENDIX A: FREE ENERGY OF THE VORTEX LATTICE STATE 



The expressi on o f the free energy F is determined so that dF/dd*{r) = and dF/ds*{r) = 0, respectively, give the 
GL equations (vM) and (2.7). Thus, the free energy is written as 



with 



+ Vo { - MP + (cs7)"^ |sp + d*Kd + s*K,sS + s*Ksdd + d*KsdS 



F . H(r)^ 



(Al) 



(A2) 



where 



i?3 = |(<^^)Mpd - \ieki>\m^ ■ q)'rf}MP + rf'{(v • c,fdY - 2d\(y ■ q)dp + 3{(v • q)4'd*]), 



(A3) 



^5 = -37(0!^) Ml V 



(A4) 



(A5) 



(A6) 



From dF/dA{Y) = 0, the expression of the current density is reproduced as follows. 



j(r)=eVxh(r)/(0o/e') 



d -——d + s 



2k2^ Uo V dA 



dA 



S + S -^r-. — d + d -^r-. — S 



dA 



1 „ai?3 

rd 



OA J 24 dA 



(A7) 



(A8) 



Equation (A8) is equivalent to Eq. (3^), if we neglect the term of d* (dRs/dA) which is the order O(MoP)- From 
d*{dF/dd*) = and s*{dF/ds*) = 0, we obtain the following relations, 



% 'MP = Vo' {d*Kd + d*Ksds) + vo^ \ d*R3 + -{fyl) (2|sPMP + s^d*^) I + Vo^'d* 



i?5 



(A9) 



Cs7??o 



% ' {-s^K^^s + s*Ksdd) + vo" [li^tM' + li^l^'^'d) (2MPkP + d^s*^ 



By substituting Eqs. (A9) and (AlO) into Eq. (Al), the free energy is reduced to 

F 



Fo 



(AlO) 



(All) 



By using the magnetic flux density B in Eq. (4.2) 
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H(r) = Ho + h(r) = B - h(r) + h(r). 



(A12) 



Thus, H(r)2 in the first term of Eq. (All) is written as 



H(r)2 =B2-h +h2. 



(A13) 



We now consider the behavior in the immediate vicinity of Hc2- The quantities d{r), s(r), A(r) and H(r) are, 
respectively, divided into the values at Hc2 (we denote them as dc, Sc, Ac and Hc2) and the deviations from them. 



d — dc + 5d, 



Ss, A = Ac + (5A, 



(A14) 



where we have relations. 



H(r)=Hc2 + {Ho-Hc2 + h(r)}, 



V X Ac = Hc2, V X (5A = Ho - Hc2 + h(r). 



By using Eq. (A14), the kernels K, K^g and Kgd are also divided into two parts, for example. 



K{A) = k(^) + ■ SA, 



where we write 



if = if(Ac), 



ax^^' dK 



dA 



dA 



(A15) 
(A16) 

(A17) 

(A18) 



While the deviations Sd/dc, 5s/ Sc and |(5A|/|Ac| are in the order of small quantity (iJo — Hc2)/Hc2, the order 
parameters dc and themselves are also small quantities in the immediate vicinity of -ffc2- The order parameters dc 
and Sc, which have been already obtained in Sec. ||, satisfy the linearized GL equations, 



-dc + K^'Uc + K^fJsc^O, 



(A19) 



ic.-f)-\sc + Ki'^^Sc + Klpc = 0. 



(A20) 



Equa tions ( A19 ) and ( A2C| ) are equivalent to Eqs. ( ^.10 ) and (2.11), respectively. By substituting Eqs. (A14) and 
( [AI^ ), Eq. (|A8|) is reduced to 



j(r)=eVxh(r)/(0o/^')-- 



d^—— dr + S,—— Sc 



dA 



dA 



^ dA ' 



dK 



sd , 



dA 



(A21) 



within t he or der 0(|doP). In this order, j(r) and h(r) in Eq. (A21) are reduced to the ones obtained in Sec. |T| [Eqs. 
( p^ , ( p^ ) and (pgj)]. 



On the other hand, by substituting Eqs. (A14) and ( A17 ) into Eqs. (A9) and ( AlO ) and by using Eqs. ( A19 ) and 
([A2C|), we obtain 



where V x h comes from the relation of Eq. (A21). From Eqs. ( A16| ) and (1.2), there is a relation 



(V X h) • 5A = (V X (5A) • h = (Ho - Hc2 + h) • h = (B - Hc2) • h - h + h^. 
Thus, by using Eqs. ( [A2^ ) and ( [A2^ ), we obtain an identity 



(A22) 



(A23) 



\do\' = 



, (B-Hc2)-h 



DAdc, Sc] 



(</>0/C2)2|^^|4 2%4|rfo| 



(A24) 
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which determines the amplitude ofthe orde r par ameter. E quati on (A24) corresponds to the "first Abrikosov identity". 
By substituting Eqs. (A13), ( A14 ) and ( A17 ) into Eq. (All), the free energy is written as 



F 



B2 



Moi'+o(idon. 



(A25) 



Equation (4J.) is obtained by the substitution of ( |A24 ) into Eq. ( A25 ). 

As for D4 de fined in Eq. (A2), the terms e xcep t for dR^ is in the order 0(7^), si nce s is in the order 0(7) as is 
seen from Eq. ( |2.26|) . Thus, D^dcSc] of Eq. (jij) is obtained by substituting Eqs. and { ^lEj ) into Eqs. (^) 

and (|A3|). 
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FIG. 1. The configuration of tlie vortex lattice and the coordinate system. The vortex centers are shown by solid circle. 
And ri and r2 are unit vectors. We set the x axis along the direction of ri . The a axis (or the 6 axis) forms an angle of 60 
from the x axis (or the y axis). 

FIG. 2. Spatial variation of the induced ,s-wavc order parameter for the case of a square lattice {(ly/dx = 0.5, C, = 0.5). 
(a) Contour plot of the amplitude, |s(r)/7S2"'' |. The region l.'dax x l-9ax is presented. The a axis and the b axis are along 
the horizontal and the vertical directions, (b) The corresponding position for the vortex of the s-wave order parameter and 
its winding number are schematically presented. The solid circle shows the vortex center of the dominant d^2_y2-wave order 
parameter, where the induced s-wave component has a vortex with winding number —1. The s-wave component has an extra 
vortex at the point marked x. The solid line shows the unit cell of the vortex lattice. 



FIG. 3. The same as Fig. 2, but for the case of an oblique lattice {ay/ax = 0.676, C = 0.5). 

FIG. 4. The equilibrium state of the vortex lattice is schematically presented for the case (a). Vortices form a shape of an 
isosceles triangle with OA = AB, where the base OA is along the a axis. 

FIG. 5. Temperature dependence of the unit cell shape of the vortex lattice. The ratio ay /ax along the curve of Hc2{T) is 
presented as a function of T/T^ for various mixing of the s-wave component, Cg = —0.2, —0.1, 0, 0.2, 0.4. A thick line shows 
the pure dj.2_j,2-wave case (cs = 0). The angle B in Fig. 4 gradually varies from 60° [ay/ax = 0.866) to 90° {ay/ax = 0.5), 
when temperature is lowered from Tc. 

FIG. 6. Spatial variation of the d^2_y2-wave order parameter. Contour plot of the amplitude, \d{r)/do\, is presented, (a) 
The case of a square lattice {ay/ax — 0.5, ( — 0.5) with an equilibrium orientation 60 = 0° for 7 — 0.2. (b) The case of a 
square lattice with an unstable orientation ^0 ~ 45° for 7 = 0.2. (c) The case of an oblique lattice (ay/ax = 0.676, ^ = 0.5) 
with an equilibrium orientation do = 0° for 7 = 0.115. The region l.Qax x 1.9ox is presented, where one of the vortex centers 
locates at the center of the figure. The a axis and the b axis are along the horizontal and the vertical directions. 



FIG. 7. Spatial variation of the current. Contour plot of the amplitude, |j(r)|/(7r€^''^^/2/t^77^), is presented, (a), (b) and 
(c) correspond to the cases of Fig. 6 (a), (b) and (c), respectively. 



FIG. 8. Spatial variation of the induced magnetic field. The contour plot of hz{r)/\hi \ is presented, (a), (b) and (c) 
correspond to the cases of Fig. 6 (a), (b) and (c), respectively. 



FIG. 9. The magnetic field distribution function P{h) as a function oi h/\hi \ for various temperatures; 7 =0 (a), 0.06 (b), 
0.115 (c), 0.15 (d) and 0.2 (e), which correspond to T/Tc =1, 0.936, 0.881, 0.848 and 0.802, respectively. Sohd lines are for the 
case of a dx2_y2-wscve superconductor, where (a) is the triangular lattice, (b) and (c) the oblique lattice, (d) and (e) the square 
lattice. There, parameters are set as C = 0.5, ^0 =0 and ay/ax =0.866 (a), 0.799 (b), 0.676 (c), 0.5 (d), 0.5 (e). Dotted lines 
are for the case of an isotropic s-wave superconductor, where (a)-(e) are the triangular lattice (^ = 0.5, ay /ax = 0.866). For 
clarity, the curves have been sifted along the vertical axis. 
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